Can a double stranded DNA be unzipped by pulling a single strand?: Phases of 

adsorbed DNA 
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We study the unzipping of a double stranded DNA (dsDNA) by applying an external force on a 
single strand while leaving the other strand free. We find that the dsDNA can be unzipped to two 
single strands if the external force exceeds a critical value. We obtain the phase diagram which is 
found to be different from the phase diagram of unzipping by pulling both the strands in opposite 
directions. In the presence of an attractive surface near DNA, the phase diagram gets modified 
drastically and shows richer surprises including a critical end point and a triple point. 

PACS numbers: 82.37.Rs, 68.35.Rh, 05.70.Jk 



I. INTRODUCTION 



DNA replication in prokaryotes gets initiated by un- 
zipping of a few base pairs at one end of the dsDNA 
and then continues till end. These processes are assisted 
by various enzymes, often by exerting force on DNA Q. 
It is now known theoretically, that a dsDNA undergoes 
an unzipping transition under the action of an external 
force if the_force_exceeds a temperature dependent crit- 
ical value 
transition 




Many studies of this unzipping 
have revealed the importance of 
ensembles namely, the fixed distance and the fixed force 
ensemble 
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For the fixed distance case, the distance 
between the strands is kept constant and the force re- 
quired to maintain the distance is allowed to fluctuate, 
whose average is the quantity of interest, whereas, the 
average of the fluctuating distance between the strands 
where an external force is acting is the quantity of inter- 
est in the fixed force ensemble. For the single molecule 
studies, the results are known to depend on the ensem- 
ble used 0, [lH- The helicases work in both the fixed 
distance and the fixed force ensembles (l3j . 

In recent years, the unzipping transition has been 
studied in detail with various extensions of the basic 
model. These include studies of models with interme- 
diate phases [f| QjJ EH EH, dependence on pulling di- 
rections fl6|. models with additional features like semi- 
flexibility ^13] , heterogeneity [1, El El H9, EH , satura- 
tion of hydrogen bonding [22j|, random forcing [23| etc. 
A similar problem of unzipping of an adsorbed polymer 
from the surface has also been studied 0, Ex- 
perimental studies use various micro-manipulation tech- 
niques [H [27], In all of these studies, the focus was 
on the breaking of the pairing by the external force. A 
more complicated situation emerges in DNA replication 
and segregation (e.g. in Bacillus subtilis, E. Coli etc.), 
where the membrane-DNA complexes play an important 
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role [29|, [30|. For example, in Bacillus subtilis, the DnaB 
protein has been shown to be a membrane-associated pro- 
tein that is involved in initiation of replication [3l|. Also, 
it is known that the interaction between the DNA and 
the membrane can form ordered domains [32| . Analo- 
gous laboratory situation would involve a substrate-DNA 
interaction. Recently atomic force microscopy has been 
used to identify the binding mechanism between the DNA 
and the DNA-binding agents by pulling one end of the 
dsDNA, which is immobilized on a gold surface in the 
presence of DNA-binding agents [33[ . 

In this paper our aim is to study a simple model of 
unzipping of a dsDNA by applying an external force on 
a single strand in the transverse direction. Can a dsDNA 
be unzipped? Within the same framework, we also study 
the behaviour of the dsDNA in the presence of an attrac- 
tive surface near it. This surface can mimic a membrane 
on which the DNA can be immobilized. Since adsorption 
of a polymer (here a dsDNA or a single strand of dsDNA) 
is a phase transition by itself (often critical), the compe- 
tition between adsorption and unzipping is expected to 
add new features to the unzipping phase diagram. 

The presence of the additional DNA substrate interac- 
tion both in-vivo and in-vitro system requires a different 
extension of the hitherto used models for the unzipping 
transition. A helical structure could pose extra prob- 
lems for surface adsorption. In order to focus on the 
competition between adsorption and pairing, we use an 
extension of the Poland-Scheraga (PS) model [34| in two 
dimensions. Previous studies of DNA unzipping showed 
that the lattice model preserves, even in two dimensions, 
the basic results of DNA unzipping including the first 
order nature of the phase transition and the existence 
of a reentrant region For the problem at hand, we 
have also done Monte Carlo simulations in 2+1 dimen- 
sions, though without considering the helical structure 
of DNA, but the strands can wind together because of 
the dynamics, and find that the force-distance isotherms 
arc qualitatively similar to the isotherms obtained in 1+1 
dimensions. 

At this point it is worthwhile to mention a few caveats 
of our approach. To achieve simplicity and tractability, 



wc ignored the helical structure of dsDNA. The main ob- 
stacle in incorporating it is the absence as yet of any 
analytically tractable models for DNA melting that ad- 
mits a helical ground state. In this situation, a PS type 
model serves as the starting point. We will see that if the 
base pairing energy of the DNA is greater than the bind- 
ing to the surface, the unzipping of dsDNA by pulling a 
single strand takes place when it is away from the sur- 
face. Since the other strand of the DNA is left free, it 
can unwind itself from the pulled strand and the helical 
structure of the DNA can be safely neglected. However, 
in the opposite limit, the helical structure of the DNA 
becomes more and more crucial as the strength of bind- 
ing to the surface is increased and its role can no longer 
be neglected. Another feature that may play an impor- 
tant role is the difference in the length per base pair of 
dsDNA and single-stranded DNA (ssDNA). Since it is 
difficult to incorporate this in a PS type lattice model 
considered in this paper, but can be done in continuous 
models, we ignore it in the present study. 

Among the new features we find is the existence of a 
critical end point (CEP) [3{|, and the triple point [IH 
in the force versus temperature phase diagram for suffi- 
ciently strong attraction with the substrate. For weaker 
attraction, certain phases may not be thermodynamically 
stable. Some of the details of the CEP have been reported 
in a shorter paper [37| . Here we give all the necessary de- 
tails and calculations and focus on the triple point. We 
use the generating function, the exact transfer matrix 
and Monte Carlo simulations to explore the equilibrium 
behaviour of this DNA-substrate system under a force on 
one strand. 

The paper is organized as follows: In Sec. [Ill we define 
our model. Section Hill is devoted to two extreme limits 
in which the problem reduces to the unzipping of ad- 
sorbed polymer from a straight and a zig zag hard-wall. 
These results are used in the subsequent sections. As 
a prelude to the substrate interaction problem in hand, 
we also need to consider the case of unzipping of a ds- 
DNA strand. In Sec. IIV1 we review the unzipping of 
a dsDNA by pulling its strands in opposite directions. 
The unzipping by pulling a single strand is studied in 
Sec. El In Sec. IVIi the unzipping of an adsorbed DNA 
by pulling a single strand is studied. The existence of 
a triple point [36| and a critical end point are shown 
here. There occurrences are not dependent on each other. 
These emerges as the relative attraction of the substrate 
is made stronger. In addition to providing some details 
on CEP this section focuses on the triple point. Finally, 
wc draw our conclusions in Sec. IVII1 
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FIG. 1: Schematic diagram of a dsDNA adsorbed on the 
surface (shaded region) in 1 + 1 dimensional square lattice. 
There is a binding energy —e^ between bases (shown by filled 
circles) of the strands of the DNA. One end of the DNA is 
always kept anchored on the surface at the origin (shown by 
star). The free strand (denoted by 1) can gain energy — e w 
for every contact with the surface (i.e., x\ = 0). An external 
force g (shown by arrow) is applied at the free end of the 
pulled strand (denoted by 2) in the transverse direction. 



allowed to cross each other but whenever they meet (i.e. 
x\ = x 2 ) there is a gain in energy — e&(e6 > 0) for every 
contact. At the diagonal (x = 0) there is an impenetra- 
ble attractive surface, with an energy — e w (e w > 0), which 
favors adsorption of the DNA. In 1 + 1 dimensions, the 
surface is a line passing through the diagonal of a square 
lattice, and only one of the strands can get adsorbed on 
it (i.e. x\ =0), since the two strands of the DNA cannot 
cross each other. One end of the DNA is always kept an- 
chored at the origin. We apply an external force <?, along 
the transverse direction (^-direction) on the free end of 
one of the strands of the DNA. The other strand is left 
free. Henceforth, the strand which is left free is called 
the "free strand" and the strand, on which the exter- 
nal force acts is called the "pulled strand" . A schematic 
diagram of the model in I + f dimensions is shown in 
Fig.[U In 2 + I dimensions, the surface is a plane passing 
through the diagonal of a cubic lattice. Unlike the I + 1 
dimensional case, both the strands of the DNA can get 
adsorbed on the surface and still satisfy the non-crossing 
constraint on the plane (y-dircction) . In this paper we 
concentrate only on 1 + 1 dimensional case. 

Let D n (x\,x 2 ) be the partition function (temperature 
dependence not shown explicitly) of a dsDNA in the fixed 
distance ensemble where nth monomers of the strands are 
at positions X\ (free strand) and x 2 (X2 > x±) (pulled 
strand) respectively from the wall. D n (x\, x 2 ) satisfies 
the recursion relation 



D n+1 (xi,x 2 ) = D n (x 1 +i,x 2 +j) 

x[l + WS Xu0 }[l+B5 X2 , Xl ], (1) 



II. MODEL 

We model the DNA by two directed self avoiding walks 
(DSAWs) on a D = 1 + 1 dimensional square lattice. The 
walks, starting from the origin, are directed along the di- 
agonal of the square (z direction). The walks are not 



where W = {e^ - 1), B = {e t)e " - 1) and (3 = 1/T 
is the inverse temperature in units of fee = 1, with the 
initial condition D (xi, x 2 ) = (e /3ew S Xlt o)(e^ 6b 5 X2t o). The 
impenetrability of the surface and the non-crossing of the 
strands demand x 2 > x\ > 0. The canonical partition 
function with an external force g at the end of the pulled 
strand is then obtained by summing over all the allowed 
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configurations of the DNA of length N on the lattice. 

Z N ((3,g)= D N ( Xl ,x 2 )e^, (2) 



X2 >X1 >0 



where e@ 9X2 is the Boltzmann weight due to the force g. 

From x\ and x%, one can define the relative coordinate 
x, and the center of mass coordinate X as 



x 2 —xi, 



X= 1 -[x 1 



■ x 2 \ 



(3) 



The relative coordinate contains all the necessary in- 
formation of the unzipping transition of a dsDNA when 
its strands are pulled in the opposite directions [1, U H @ ■ 
For the single strand pulling case we need both the rela- 
tive and the center of mass coordinate to track the chains 
individually. 



A. Quantities of Interest 



The following physical quantities are of interest: 

1. The average distances of the end monomers of both 
strands from the surface, which arc defined as 



E 

X-2>X\ >0 



Xi D N { Xl ,x 2 ) 
Zn(P,9) 



(4) 



where (• • • ) denotes the thermal averaging and sub- 
script i=l and i = 2 stand for the free and the 
pulled strand respectively. The summation is taken 
over all the allowed configurations. The positions 
of end monomers give us the information needed to 
characterize the phase of the system. 

2. The response to the force, i.e. the isothermal exten- 
sibility, which can be expressed, for both strands, 
in terms of fluctuations of the position of end 
monomers 



d(xi 



dg 



1 



(i = 1 or 2). 



(5) 



Equations ([4J and ([5]) are useful to obtain the above phys- 
ical quantities numerically. 



mass frame. The problem then reduces to the unzipping 
of an adsorbed polymer on an straight impenetrable sur- 
face (hard- wall). However, for case (ii), the free strand 
remains adsorbed on the surface and itself acts as a zig 
zag hard-wall for the pulled strand. In this section, we 
obtain the phase diagram for both the cases. 



A. Straight hard- wall 

In the fixed distance ensemble, the partition function 
w n (x) of a dsDNA in the center of mass frame satisfies 
the recursion relation 

w n+1 (X) = [w n (X + 1) + w n (X - 1)] [1 + WS x ,o] . (6) 

Note that while writing Eq. (O from Eq. (JTJ) we have 
absorbed the Boltzmann factor e@ eh in w n (X). The par- 
tition function in the presence of an applied external force 
g is then given by 



X>0 



w N (X) e 



(I) 



This problem has been studied in recent years because 
of its similarity with the DNA unzipping [24|, [2{|. It 
is known that the polymer unzips from the surface if the 
force exceeds a critical value. Below this critical force, the 
polymer remains adsorbed on the surface while above it 
the polymer is in the desorbed phase. The phase bound- 
ary separating the two phases is given by 



g^(T) = Un [e" e 



T 



(8) 



and is shown in Fig. [2] by a solid line. The triangles on 
the line represent the phase boundary obtained numeri- 
cally by using the exact transfer matrix. The method is 
introduced below for the zig zag surface. Here we have 
taken e w = 1. The critical force decreases monotonically 
with the increase of temperature and becomes zero at 



T C W =e w /ln2, 



(9) 



where the polymer desorbs from the surface because of 
the thermal fluctuations. 



B. Zig zag surface 



III. UNZIPPING AN ADSORBED POLYMER 

Let us first study the two extreme limits: (i) £b — > oo 
and e w is finite (and we take e w = 1), and (ii) e w — > oo 
and eb is finite (eb = 1). For the former case, the two 
strands of the DNA always stay together for the entire 
range of temperature in which we are interested. In such 
a situation, an external force g on the pulled strand also 
pulls the free strand. Therefore, the DNA can be equiv- 
alently represented by a single polymer in the center of 



The recursion relation for the partition function b n (x2) 
of the pulled strand, in the fixed distance ensemble, can 
be written as 

b„+i(x 2 ) = [b„(x 2 + 1) + b n (x 2 - 1)] 

x [1 + B6 Xat0 ] [1 + B6 X2tl ] . (10) 

Note that the above recursion relation differs from 
Eq. © in the extra factor due to the binding of the 
odd monomers to the wall. In Eq. only the even 
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FIG. 2: g vs T phase diagram for the unzipping of an ad- 
sorbed polymer from an attractive impenetrable surface. The 
points are obtained by using the exact transfer matrix and the 
solid curve is an analytical result (see Eq. (jHJ). The triangles 
and squares represent the phase boundary between the ad- 
sorbed (Ads) and the desorbed (Des) phases for the straight 
and zig zag hard-walls respectively. The thermal desorption 
temperatures for the straight and the zig zag surfaces are 



shown respectively by T c 
shown by a circle. 
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monomers bind to the surface. This modification is 
enough to make a difference in the phase diagram. This 
extra energy affects the intermediate temperature be- 
haviour of the polymer as it becomes energetically costly 
to create bubbles on the polymer adsorbed on a zig 
zag hard-wall. This can be understood as follows: The 
ground state for the polymer adsorbed on a straight hard- 
wall contains inherent bubbles of length 1 = 2 because 
the geometry of the problem allows only N/2 monomers 
on the wall. In contrast, for the zig zag hard- wall there 
are no such inherent bubbles as all the TV monomers are 
adsorbed on the wall. To create a bubble of length £ = 2, 
one of the monomers has to come out from the wall which 
costs energy. Furthermore, getting larger bubbles are 
easier for the polymer adsorbed on the straight hard- 
wall than the zig zag case. For example, flipping of one 
monomer from the wall can create a bubble of length 
I = 4 for the normal hard-wall, while three consecutive 
monomers have to come out from the wall for the zig zag 
case to create a bubble of the same size. 

Under the influence of a fixed pulling force g at the 
free end, the partition function is obtained by 

Z N (P,g)=J2b N (x 2 )e' 3 ^. (11) 

The partition function of the chain length N is ob- 
tained numerically, at a given temperature, by iterating 
Eq. pop and, then, for various force g by using Eq. (fTTj) . 
This is known as the exact transfer matrix technique. 
The distance of the end monomer of the pulled strand 
from the surface, which is the quantity of interest, is ob- 
tained by using 

(X2) =77UE*N ete ' ( 12 ) 



FIG. 3: Data collapse of the average distance (22) of the last 
monomer from the zig zag surface (a) as a function of g at 
T = 1.0 for TV = 1000, 2000 and 3000. The critical exponents 
are d = 0.99 ± 0.01 and cf> = 1.0 ± 0.01 and the critical force 
gi z) = 1.059 ± 0.001. (b) as a function of T for N = 2000, 
3000, 4000 and 5000. In this case the critical exponents are 
d t = 0.53 ± 0.03, 4> t = 0.4 ± 0.03 and T c (z) = 4.81 ± 0.16. 



In Fig. EJa), we have shown the data collapse of g vs 
(X2) isotherms at T = 1.0 for the chain of lengths N = 
1000, 2000 and 3000. We have used the scaling form 

(x 2 ) = N d y ((«? - g^)N*) , (13) 

(z) 

where d and 4> are the critical exponents and <7c is the 
critical force. By using the Bhattacharjee-Seno proce- 
dure [Hi, we obtained g [ c z) = 1.059±0.001, d = 0.99±0.01 
and 4> = 1-0 ± 0.01. These exponents are same as the 
exponents obtained for the unzipping from the straight 
hard- wall |24| and the DNA unzipping problem [J@|. We 
use same procedure at various temperatures to obtain the 
phase diagram. 

The phase diagram of unzipping from the zig zag sur- 
face is shown in Fig.[^]by points for ej, = 1. If we compare 
this with the phase boundary for the straight surface, we 
see that at very low temperatures (T = 0), the force 
required to unzip the polymer from the zig zag surface 
is twice that from the normal surface since one has to 
overcome an additional binding at a; 2 = 1. This implies 
the low temperature behavior of the polymer is same for 
both the walls i.e. g on e ( e = e w for the straight surface 
whereas e = 2^ for the zig zag case) . The polymer starts 
behaving differently as temperature is increased. As al- 
ready mentioned, it is easy to create a bubble on the 
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FIG. 4: g vs T phase diagram for the DNA unzipping. The 
dashed line represents the phase boundary for the unzipping 
by pulling strands in opposite directions by an external force 
force g. The solid line is exact phase boundary (Eq. |)24p 1 
and the triangles are from the numerics for the pulling of a 
single strand for which the transition takes place only above 
a minimum temperature T u . 



polymer adsorbed on the straight wall, whereas more en- 
ergy is needed to create a bubble of the same size on the 
polymer adsorbed on the zig zag surface. For the straight 
surface, the critical force needed to unzip the polymer 
from the wall decreases monotonically with the temper- 
ature. In contrast, for the zig zag surface, the critical 
force increases at intermediate temperatures, reaches to 
a maximum value and then decreases to zero, thus show- 
ing reentrance as shown in the inset of Fig. [3] To obtain 
the desorption temperature for the zig zag case, we resort 
again to the finite size scaling (but with temperature as 
a variable) of the form 

(x) = N dt y ((T - TW)N+A , (14) 

with Tc Z ^ as the critical desorption temperature and d t 
and fa are the critical exponents. The data collapse, for 
chain lengths N = 2000, 3000, 4000 and 5000, obtained 
for d t = 0.53±0.03, fa = 0.4±0.03 and T c (z) = 4.81±0.16 
is shown in Fig. [3{b), and the desorption temperature 



Tc is shown in Fig. [5] by a solid circle. 



IV. UNZIPPING DNA BY PULLING STRANDS 
IN OPPOSITE DIRECTIONS 



Before considering the unzipping by pulling a single 
strand, let us first concentrate on the unzipping of a ds- 
DNA by pulling its strands in opposite direction by an 
external force g. This problem has received a lot of at- 
tention in recent years due to its resemblance with the 
way the DNA unzipping experiments are done [2ol l28j . 

In the fixed distance ensemble, the partition function 
of the dsDNA can be obtained by assigning e w = in 
Eq. {!]). The recursion relation, in the relative coordi- 



nate, reads as [6| 

d n+1 (x) = [d n (x + 1) + 2d n {x) + d n (x - 1)] [1 + B6 x ,o] , 

(15) 

with the initial condition do(x) = e^ tb 8 x .Q. 

In the fixed force ensemble, the partition function for 
the DNA of length N, is then obtained by 



Z N (f3,g) = y^djy(x) 



Jig* 



(16) 



The phases of the DNA and the transition come from 
the singularities of the of the generating function 



g(z,f3,g)=f^z N Z N (l3 > g). 



N=0 



The singularities are 0, 0] 



Zl 



-~2 



09, cb) = V 7 ! - er^ - 1 + e 



0£h 



and 



1 



2 + 2 cosh (3g 



(17) 



(18a) 



(18b) 



(18c) 



The phase of the DNA is given by the singularity which 
is closest to the origin and the phase transition takes 
place when the two singularities cross each other. For low 
force, Z2Q3, £b) is closest to the origin and the DNA is in 
the zipped phase (double-stranded) , while for high force, 
z-j,{J5 1 g) becomes closest and the DNA is in the unzipped 
phase (two single strands). The force-temperature phase 
boundary between the two phases is obtained by equating 
the two singularities, which gives 0, Q 



T 

9c(T,e h ) = —cosh 
T 



2z 2 (P 1 e h ) 



- 1 



= -^^\(z 2 (f3,e h )) , 



(19) 



where X(z) = (1 - 2z - y/l - 4z)/(2z). The thermal 
denaturation (melting), coming from z\ = z 2 (P, £b), is at 



e b /ln(4/3). 



(20) 



The phase boundary separating the zipped phase (Zp) 
from the unzipped phase (Uz) is shown by the dashed 
line in Fig. [4j 



V. UNZIPPING DNA BY PULLING A SINGLE 
STRAND 

In previous sections, we have seen that the finite size 
scaling of the force-distance isotherms (or the extensi- 
bility), obtained by using the exact transfer matrix for 
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various chain lengths, can be used to obtain the phase di- 
agrams numerically. We use Eqs. |T]) and @ with e w = 
to obtain the partition function of a dsDNA of length N. 
The average distances of the end monomers of both the 
strands from the surface is obtained by using Eq. ([4]). 

I 



with the initial condition do(x) = e^ eb S x .Q. The above 
recursion relation can be analyzed analytically and the 
boundary separating the zipped and the unzipped phases 
can be obtained exactly. The details of the calculation 
are given in Appendix A. 



A. Isotherms and Extensibility 

The force-distance isotherms at T = 0.5 and 1.5 are 
shown in Fig. OJa). Due to the entropic repulsion, the 
dsDNA stays at a distance of y/N from the surface, even 
for g = 0, to maximize its entropy (see Fig. [SJb)). When 
the pulling force, <?, is small, the binding energy, e^, wins 




0.02 0.04 0.06 0.08 0.1 
(b) g 

FIG. 5: (a) g vs (xi, 2 )/N isotherms at T = 0.5 and 1.5 for 
e w = 0. The squares and circles refer to the two different 
strands of the dsDNA. The critical force above which the 
dsDNA unzips is shown by g u - In (b) same isotherms are 
shown near g = in a different scale (xi^) /N 1 ^ 2 . The broken 
lines are a guide to the eyes. 



Alternatively, using the fact that the force is only on 
one of the strands and the surface plays no role in the 
phase boundary, the partition function can be obtained, 
in the relative coordinates, by the following recursion re- 
lation (in a mixed fixed-distance-force ensemble) 



(21) 

I 

over the entropy, which the free strand can gain if it 
separates itself from the pulled strand, and the dsDNA 
as a whole gets stretched in the direction of the force. The 
average distances of end points of the strands from the 
surface, (2:1,2) /N, remain the same (for both the strands) 
and increase linearly with g. The slope, however, depends 
on the temperature; it is larger at low temperatures and 
smaller at high temperatures. As g is increased further, 
the dsDNA gets completely stretched at T = 0.5. But 
for T = 1.5, before the dsDNA can get fully stretched, 
a critical force </ u is reached, and the free strand of the 
DNA gets unzipped from the pulled strand to increase its 
entropy. We call this as "transition Uz" . This transition 
can be studied in the relative coordinate, x, defined in 
Eq. ([3]). Below g u , the strands stay together {(x)/N — > 
as ./V — > 00), and above it, they are maximally separated 
((x)~N). 

The isothermal extensibility, obtained by using Eq. I|5*|). 
is plotted in Fig. E^a) as a function of g for various chain 
lengths at T = 1.5. The critical force g u can be located 
by using the finite size scaling of the form 

X = N d G((g-g u ) N+), (22) 

with d and <j) as critical exponents. By using the 
Bhattacharjee-Seno procedure, we obtained d = 2.02 ± 
0.01, = 1.01±0.01 with g u = 2.692±0.001. These expo- 
nents are same as the exponents for the DNA unzip pin g 
by pulling strands in the opposite directions [J, la. Il0j|. 
This indicates that, similar to the later case, the unzip- 
ping by pulling a single strand is also a first order phase 
transition with 

X/N ~ \g-g u T 1 . (23) 

The data collapse is shown in Fig. [6jb). The finite size 
scaling of extensibility, as described above, can be used 
at various temperatures to obtain the phase diagram of 
unzipping. 



B. Phase Diagram 

The boundary separating the zipped (Zp) and the un- 
zipped (Uz) phases is given by (see Appendix A for de- 



dn+i(x) = [d n (x + 1) + d n {x) eP a + d n (x) + d n (x ~ 1) e'^] [1 + BS X>0 ] , 
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FIG. 6: (a) Isothermal extensibility \u vs g at T = 1.5 for 
the DNA of lengths N = 2000, 2500 and 3000. (b) The data 
collapse of the extensibility. For all plots 6b = 1. 



tails) 



fc B T, /2e" /3eb 
In 



1 - 2e-^ 



(24) 



The same phase boundary has also been obtained by 
Marenduzzo et al. \^ in a different context of melting of 
a stretched DNA. The phase boundary is shown in Fig. 0] 
by a solid line which matches excellently with the results 
(triangles) obtained numerically from the exact trans- 
fer matrix. The plot shows that in contrast to pulling 
both the strands in opposite directions, where, below the 
melting temperature, the DNA can be unzipped at all 
temperatures including T = 0, the DNA can only be 
unzipped above a certain temperature T u (T u = en/ In 2) 
for the single strand pulling. Below T u , the DNA remains 
in the zipped phase for any value of force g. The melt- 
ing temperature of the dsDNA remains the same, i.e. 
T m = eb/ln(4/3), because in our model the force acts 
in the transverse direction and does not overstretch the 
DNA. Recently it was found that a longitudinal stretch- 
ing force on one of the strand destabilize the DNA which 
results in a reduced melting temperature SH IH, EI| • 
The origin of the unzipping transition by pulling a sin- 
gle strand is different from the unzipping transition by 
pulling both the strands. For the latter case, the tran- 
sition sets in due to the competition between the base 
pairing energy et,, which binds the complementary bases 
(or monomers) of two strands, and the orientation of 
the individual links connecting the monomers. In con- 
trast, the interplay between the binding energy et,, and 



the entropy, which the free strand can gain if it is in the 
unzipped phase, is responsible for the single chain case. 
When T < T u , the binding energy wins over the entropy 
and the DNA remains in the zipped phase for any value 
of g. For large g, the DNA takes a fully stretched configu- 
ration and bubbles are not possible. This is analogous to 
the Y model studied in Ref. @ . In this model, the ther- 
mal melting of dsDNA takes place at T u = et,/ In 2 and 
the transition is of first order Q. For T < T u , the DNA 
always remains in the zipped phase. For T > T u , the 
free strand has all the favourable conditions to increase 
its entropy. Therefore, as T is increased, the free strand 
can get separated from the pulled strand well below the 
fully stretched configuration and the critical force g u falls 
rapidly with increasing the temperature, becoming zero 
at T m . 

There are other factors that also contribute in the un- 
zipping of dsDNA by pulling a single strand but neglected 
in our lattice model. An important one is the fact that 
the length per base pair of the dsDNA and the ssDNA 
arc different (0.34nm and around 0.5mn respectively). If 
the applied force is low, the dsDNA has larger extension 
than the ssDNA and the dsDNA remains stable. How- 
ever, for sufficiently strong force, the extension of ssDNA 
becomes more than the dsDNA and the force destabilize 
the dsDNA and favors its unzipping [40]. It is difficult 
to incorporate this in a lattice model like ours, but we 
believe that it can be done in a continuous description of 
the model. It would be interesting to study the combined 
effect of both mechanisms in DNA unzipping by pulling 
a single strand as this could bring down T u . 



VI. UNZIPPING AN ADSORBED DNA 

Let us now consider the complete model of an attrac- 
tive surface (i.e., e w > 0) near a dsDNA. Here the free 
strand of the DNA experiences two different energies of 
opposite tendencies even at T — 0. The energy e w tries 
to keep the free strand adsorbed on the surface while 
the binding between the strands, eb, tries to keep it with 
the pulled strand. For T > 0, the entropy also plays a 
role and the competition among them makes the phase 
diagram very rich. We establish the possibility of four 
distinct phases (I) Z a : zipped DNA adsorbed on the sur- 
face, (II) Zd: zipped DNA desorbed from the surface, 
(III) U ac i: unzipped DNA with the free strand adsorbed 
on the surface, and (IV) Udd: unzipped DNA with both 
the strands desorbed from the surface. 



A. Isotherms and Extensibility 

In Fig. we have shown the force-distance isotherms 
for e w = 1.0, 1.8 and 2.0 at two different temperatures 
T = 0.5 and 1.5 for the chain of length N = 1000. The 
values of e w and T are selected to display the typical 
characteristics of the phase diagrams. The phase dia- 
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FIG. 7: The (xi t 2)/N vs g isotherms at temperatures, T — 0.5 
and 1.5 for All isotherms are for N — 1000 with 

6b = 1. The transitions Uz and Sz take place at critical 
forces g u and g a respectively (see text). These are shown by 
arrows in the plot, (a) For e w = 1, (b) For e w = 1.8. We 
have also shown the averages obtained by using Monte Carlo 
simulations. The squares and circles refer to the two different 
strands of the dsDNA. The upper and the lower triangles are 
the estimates given by the multiple histogram technique at 
various g. (c) For e w = 2. The broken lines are a guide to the 
eyes. 



fully stretched configuration for T — 0.5. However, for 
T = 1.5, the free strand gets separated from the pulled 
strand (transition Uz) at g u and stays at a distance of 
\/~N from the surface (see Fig.JTJa)). 

The isotherms for e w = 1.8 is shown in Fig. [7£b). At 
T = 0.5, the isotherm is similar to that of e w = 1, though 
with a higher critical force g s . It is different for T = 1.5, 
where, on increasing the force <?, the transition Sz takes 
place at g s . On increasing g further, we have an another 
transition Uz at g u , where the free strand separates itself 
from the pulled end and gets adsorbed on the surface. We 
have obtained the isotherm at T = 1.5 by two different 
methods. The smaller symbols with broken lines between 
them are from the exact transfer matrix whereas, the 
bigger symbols (squares and circles representing the two 
strands of the DNA) arc obtained by performing Monte 
Carlo simulations, at their respective g values. We also 
collect histograms, hk{E,x\,X2) (E is the total binding 
energy, x\ and X2 are respectively the distances, from 
the surface, for the free and the pulled strand, and the 
subscript k, stands for the simulation performed at force 
<?fc) at each simulation. These histograms arc then com- 
bined to estimate (0:1,2) for a r ang e of forces by using the 
multiple histogram technique 44(. The estimates, so ob- 
tained, are shown by the upper and the lower triangles for 
the free and the pulled strand respectively in Fig. [7jb) . 
The isotherm obtained from both the methods agree ex- 
cellently The details of the Monte Carlo simulation is 
discussed in Appendix B. Our Monte Carlo simulations 
in 2 + 1 dimensions also shows similar results ■ 

The isotherms for e w = 2 are shown in Fig. [7]Jc). The 
plot shows that the only possible transition is the un- 
zipping of the dsDNA to two single strands (i.e. the 
transition Uz) at </ u , because the free strand minimizes 
its energy by staying adsorbed on the surface. 

The critical force <7 S at which the transition Sz takes 
place is obtained by the finite size scaling of isothermal 
extensibility as given by Eq. ([22"]) . A good data collapse is 
obtained for the exponents d s = 1.95±0.05 and <j> B = 1.0± 
0.01 showing that, as for transition Uz, the transition 
Sz is also first order. We use the finite size scaling of 
extensibility at various temperatures to obtain the phase 
diagrams. 



grams arc obtained by repeated use of finite size scaling 
of the response function. 

For e w = 1.0, the isotherms are shown in Fig.[7|Ja). In 
the absence of a force g, the ground state is an adsorbed 
DNA on the surface. As g is increased, we see that there 
is a critical force, g s , at which the DNA gets unzipped 
from the surface but remains double stranded. We call 
this as "transition Sz" . This involves the center of mass 
coordinate, X, given by Eq. (J3]). Below g s , (X)/N = 
as N — > 00, but it takes a finite (non-zero) value above 
g s . When g is increased further, the dsDNA stretches 
more and more in the direction of the force, and takes a 



B. Phase Diagram 

The discussions in the previous subsection reveal that 
there are four possible phases for the problem at hand. 
These are: 

1. Phase I — The zipped DNA (i.e., dsDNA) ad- 
sorbed on the surface. This phase is character- 
ized by (2:1,2) /AT -> as JV -> 00. This phase 
is called Z a . 

2. Phase II — The dsDNA desorbed from the surface. 
In this phase we have (xi^/N — > A (0 < A < 1) 
as A^ — > 00. This phase is called Z<\- 
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3. Phase III — In this phase, the DNA is unzipped 
(i.e., the strands stay away from each other). The 
pulled strand is stretched in the direction of the 
force and the free strand is adsorbed on the surface. 
This phase is characterized by (x\)/N — > but 
(x2)/N — > 1 as N — ► oo. This phase is called U a d- 

4. Phase IV — In this phase both the strands stay 
away from the surface as well as from each other. 
The pulled strand follows the pulling force as in 
phase III, and the free strand stays away from the 
surface to maximize its entropy. This phase is char- 
acterized by {x 2 )/N — > 1 and (xi)/yN — > 1 as 
N — > oo. This phase is called Udd- 

Depending on the relative strength of the binding en- 
ergy, e w , and the pairing energy eb, wc can either have 
all the four phases listed above, or a few of them, in 
the phase diagram. We have four parameters: g, T, e w 
and eb- Out of them we can construct three dimension- 
less quantities g — > g/e\>, T — > fceT/eb and e w — > e w /eb- 
Therefore, without loss of generality, we chose eb = 1. 
The phase of the DNA, for the given set of parameters, 
can be read from a 3-dimensional g-T-e w surface. Since 
it is difficult to show a 3-dimensional plot, wc show the 
cross-sections (g-T plane) of the above surface for various 
e w . These are shown in Fig. [H] 

Before discussing the phase diagrams in detail, let us 
do a zero temperature (T = 0) analysis of the problem, 
keeping constant. The energies of the three phases, 
namely Z a , Zd, and U a d are respectively given by 



and 



E Za =-N(e w /2+l), 



E Zd =-N(g+l), 



Ev ad =-N(e w /2 + g). 



(25a) 



(25b) 



(25c) 



For 1 < e w < 2, the phase U a d is always unstable (i.e. 
it has higher energy). The transition from phase Z a to 
phase Zd occurs at g = e w /2. But for e w > 2, Zd is not 
possible. At e w = 2, there is a degeneracy for Zd and 
U a d, which occurs at g = e w /2. Therefore at this point 
all the three phases coexist and it is a triple point. 



1, < e w < 1 and 1 < e w < 2 

When < e w < 1, the phase diagram contains three 
phases, namely Z a , Zd and Udd- But, for 1 < e w < 2, a 
new phase, U a d, also appears in the phase diagram. The 
phase boundary separating Z a from Zd which decreasing 
monotonically for e w < 1 now shows reentrance at in- 
termediate temperatures. Apart from this feature, there 
is a region in the phase diagram which involves three 
phases, namely Zd, U a d and U a d- The transition from 
phase Zd to U a d and Udd are of first order, whereas, the 




FIG. 8: g vs T phase diagrams for e w = 2. The points are 
obtained by using the exact transfer matrix and T w , and T m 
(shown by arrows) represent respectively the temperature at 
which the dsDNA desorbs from the surface, and the melting 
temperature of the DNA. The triple point and the critical 
end point are shown by a diamond and a square and repre- 
sented respectively by C p and C c . The thick dashed, solid and 
dotted lines are respectively Eqs. l|8]l. (|21[) and (|26p and are 
approximation to the phase boundaries represented by circles, 
up triangles and down triangles. 



transition from phase U a d to phase Udd is a second order. 
The second order phase boundary terminates on the first 
order line at a critical end point. These are discussed in 
Ref. [13. 



2. 



— 2: triple point and critical end point 



The phase diagram for e w = 2 is shown in Fig. [8] It 
contains all the four phases, namely Z a , Zd, U a d, and 
Udd, same as 1 < e w < 2 case. The phase boundary 
separating phase Z a from phase Zd is shown by circles. 
For comparison, we have also shown, by the dashed line, 
the phase boundary (see Eq. ([8]) ) for the unzipping of 
an adsorbed polymer from the straight hard- wall. At 
low T, the two strands of the DNA, which are bound 
to each other, behave like a single chain, and therefore, 
the low-T phase boundary for both the problems match 
with g(T = 0) = e w /2. As T is increased, bubbles form 
on the DNA. These bubbles create an entropic hindrance 
on the free strand to stay away from the surface. As a 
result desorption of dsDNA from the surface occur at a 
temperature higher than the temperature needed to des- 
orb a single polymer i.e. T w > e w /ln2. The boundary 
separating phase Zd with phase Udd, which is shown by 
a solid line, is again given by Eq. (f2"4"| . i.e., it is same 
as the phase boundary between the zipped and the un- 
zipped phases for e w = case. Therefore, in this region, 
the surface plays no role in determining the phase bound- 
ary between the zipped and the unzipped phases. The 
transition from phase Zd to U a d, which is well approxi- 
mated by 



9{T) 



knT 



In 



2VF^ 



-0e b 



-0e b 



-2,3e„ 



(26) 
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for the e < 2 case (sec Appendix A and Ref. (37}) is 
shown by the dotted lines. This curve deviates markedly 
with the data obtained by the exact transfer matrix cal- 
culations at intermediate temperatures, but still matches 
in both the high and the low temperatures. This devi- 
ation is because we can no longer neglect the effect of 
surface on the pulled strand mediated via the free strand 
and the approximation on which the Eq. (|26[) is based 
breaks down. At low temperature side, it meets with 
the dashed line at T = 0, showing that the three phases 
coexist there as expected from zero temperature analysis 
discussed above. This is the triple point (T p = 0) which is 
shown by C p (filled diamond) in the phase diagram. The 
numerics, however, cannot resolve the two phase bound- 
aries at low temperatures. On the high temperature side, 
the dotted curve meets with the solid line with the same 
slope but with different curvature at the CEP which is 
shown by C e (filled square). At this point the boundary 
separating phases U ac i and Udd gets terminated. 



but takes place up to the denaturation temperature. For 
both the cases, the transition is first order. 

On introducing the binding energy, e w , at the surface, 
the phase diagram becomes very rich, with a total of four 
possible phases. Depending upon the relative strength of 
e w and et>, we can either have all the four phases, or a 
few of them, in the phase diagram. We find that for a 
wide range of e w , there is a critical end point in the phase 
diagram, where a line of second order (critical line) ter- 
minates on a first order phase boundary. Furthermore, 
for a narrow range of e w , we have a triple point in the 
phase diagram. As e w — > oo, the problem reduces to the 
unzipping of a polymer, which is adsorbed on a zig zag 
hard-wall. It seems that the unzipping of an adsorbed 
dsDNA by pulling a single strand can be a potential can- 
didate to explore the critical end point in single molecule 
experiments. 

APPENDIX A: EXACT PHASE BOUNDARIES 



> 2 



When e w is slightly greater than 2, the triple point is 
at a finite temperature (T p > 0). However, it is difficult 
to get the exact location of the triple point because of the 
difficulty in resolving the phase boundaries at low tem- 
peratures. On increasing e w further, the region represent- 
ing phase Zd shrinks rapidly and both the triple point and 
the CEP shift towards higher temperature and disappear 
independently from the phase diagram. When e w = oo, 
the free strand, which remains adsorbed on the surface 
at all temperatures, acts like a zig zag hard-wall studied 
in Sec. IIIIl In this case, only the phase Z a and the phase 
U a d survive in the phase diagram as shown in Fig. [5] 
This picture is true only in D = 1 + 1 dimensions. In 
D = 2 + 1 dimensions the winding of chains will produce 
dynamic hindrance to unzipping unless torque releasing 
mechanism is provided at the anchored point. Same will 
be true for a real DNA, where its helical structure will 
prevent its unzipping though topoisomerases may help in 
releasing the extra strain. There is also the possibility of 
dynamic hindrance leading to noncquilibrium long lived 
states. These are beyond the scope of this work. 



VII. CONCLUSIONS 

In this paper, we have studied the unzipping of a ds- 
DNA by pulling a single strand. We find that a dsDNA 
can be unzipped, even by pulling a single strand, if the 
pulling force exceeds the critical value. The origin of 
this transition is different from the unzipping of DNA 
by pulling its strands in opposite directions. In contrast 
to the latter case, where the unzipping is possible at all 
temperatures below the denaturation temperature of the 
DNA, the unzipping only starts above some minimum 
temperature (T u = £b/ln2) for the single strand pulling 



In this appendix we give details of the analytical cal- 
culation for obtaining the phase boundaries by using the 
recursion relation given by Eq. (f2"Tj) . We define the gen- 
erating function of the partition function d n (x) by 



G(z,x) = ^2 zn d n (x), 

n=0 

which can be taken of the form (ansatz) 
G(z,x) = \ x {z)A(z), 



(A-1) 



(A-2) 



where X(z) and A{z) need to be determined. Using this 
ansatz in Eq. (|2"Tj) (z dependence of A and A suppressed) 
we get 



.4 

z 



{e 0g + e- 0g (l + X)}A + 



and 

A* 
z 

from which one obtains 

l r 



e 0eb for x = 0, 

(A-3a) 



-- = A"" 1 (1 + A) [e P9 + Xe-? 9 ] for x > 0, (A-3b) 



A 



e J3a _ ( e Wa 



Dz 



+ \J {(e 2 03 + I) z - ePs}- 4z 2 e 2 P9 , (A-4a) 



and 



A = 



1 - z [eJ 3( J + e- B s (1 + A)] eP £ * ' 
Now from Eqs. (|A-1|) and (|A-2|) . we have 



G{z)=Y J G(z,x)=Y, XXA =i 



A' 



(A-4b) 



(A-5) 
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which has a singularity at A = 1. Substituting A = 1 
in Eqs. (|A-4a[) and (|A-4b[) . we get singularities of the 
unzipped and the zipped phases respectively as 



1 



24 



and 



"-2 



4 cosh j3g 
2 J \2coshf3g 

e -f3e b 



(A-6a) 



(A-6b) 



The first factor in Eq. (|A-6a[) is the contribution of the 
random walk of the free strand when the DNA is in the 
unzipped phase and the second factor is the contribution 
of the force on the pulled strand. The phase boundary 
separating the zipped and the unzipped phases is given 
by Z2 = 24, which reads as 



k B T (2e-f>*> 
9{T) = -r- In 



1 - 2e~l 3 ^ 



(A-7) 



The zero force limit (i.e. g — » 0) of Eq. (|A-7|) gives the 
thermal melting temperature of the dsDNA. This gives 



T = 



In 4/3' 



(A-8) 



In the opposite limit (i.e. 
mum temperature 



00), wc obtain the mini- 



rp _ £ b 

u_ In 2' 



(A-9) 



above which the dsDNA can be unzipped by pulling a 
single strand. 

By using the same analysis given above, we can also 
approximate the phase boundaries between the zipped 
and the unzipped phases in the presence of an attrac- 
tive surface with energy — e w near DNA. By maintaining 
the structure of the singularity 2:4, the singularity which 
contribute to the partition function of phase U a( j in the 
thermodynamic limit can be approximated as 



Z3 



=,-2/3e v 



1 



2 cosh j3g 



(A-10) 



In Eq. (|A-10[) . the first factor is the contribution of the 
adsorbed free strand on the surface and the second factor 
is same as in Eq. (|A-6a[l . The phase boundary separating 
phase Zd with phase U a d is given by 



In 



2Ve^ 



-2/3e„ 



-0e b _ ^/ e -/3e» 



(A-ll) 



which is obtained by equating 2:4 with z%. This phase 
boundary is shown by dotted lines in Fig. [5] The ap- 
proximation done in Eq. (|A-10|) is valid only for smaller 
values of e w and breaks as e w becomes comparable to £b 



as evident from Fig. E^b) and (c). However, the g — > 00 
limit of Eq. (jA-llj) gives the T u and e w dependence as 



e w = k B T u In 



1- v /l-4exp(-2e b /fc B T u ) 



, (A-12) 



from which one obtains T u — > eb/m2 as e w — > 1 and 
T n as e w -> 2. 



APPENDIX B: DETAILS OF MOTE CARLO 
SIMULATIONS 

In this appendix we give the details of Monte 
Carlo simulations used in obtaining the force-distance 
isotherms in both 1 + 1 and 2 + 1 dimensions. We model 
the bases of the DNA by beads and the two adjacent 
beads on a strand, are connected by a rigid rod of unit 
length which stays on edges of the square or the cubic 
lattice in 1 + 1 and 2 + 1 dimensions respectively. If the 
beads of the two strands are unit distance apart, there is 
a binding between them. At the diagonal of the square 
and the cube there is an impenetrable surface. One end 
of both the chains are anchored at the origin and a pulling 
force g is applied on the bead at the free end of the pulled 
chain. We consider a single bead flip dynamics. In 1 + 1 
dimensions, the jth bead in the interior of the strand, 
e.g. the free strand, located at a distance x\. from the 

surface (state fx), is flipped to X-y . X 1 , ± 2 (state v), 

provided x\ j+1 — x'^. = x[. — x\ j _ 1 = 1. This constraint 
ensures that the chain does not break while doing the 
dynamics. The move is accepted (or rejected) according 
to the rule 



P{fi 



if E t 
otherwise 



> 



(B-l) 



The same thing has to be repeated for the pulled strand 
also. For the end bead we have to take care of the energy 
contribution g(x' 2N — X2 N ) due to the force g. In 2 + 1 
dimensions, the bead can be flipped in three positions 
without breaking the chain. Apart from the connectiv- 
ity constraint, the excluded volume effects, between the 
strands and between the free strand and the surface, have 
to be taken care of. This means the following moves have 
to be rejected while doing the dynamics (i) the free chain 
crossing the pulled chain, (ii) the free chain crossing the 
surface on which it is adsorbed, and (iii) the pulled chain 
crossing the first chain. It is possible for the DNA to 
reach any state from any other state using the above 
moves. The 2N such flips, N for each strand, constitutes 
one Monte Carlo step (MCS) per bead. At each value 
of the force, where the simulation is performed, we allow 
the system to equilibrate by repeating the above proce- 
dure for 10 6 MCS and start measurements only after it. 
Between any two measurements wc run the procedure, 
without doing measurements, for 10 3 MCS to avoid cor- 
relations between two measurements. 
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We have performed K — 7 simulations at j3 = 2/3 
at various g as indicated in Fig. EJb). The expectation 
value of (2:1,2) is obtained by averaging over 10 6 MCS. 
Other than the averaged values, we have also obtained 
histograms, hk(E,x\,x 2 )(k = 1, . . . K) of the energy E 
(binding from the surface plus base pairing but not the 
contribution from the force), the positions of the end 



monomers for the free and the pulled strands represented 
by x\ and X2 respectively. The fcth histogram is collected 
at force gk with nu number of states. By using the multi- 
ple histogram technique [44| . the average distance (from 
the surface) of the last monomers of both strands, (2:1,2) 
at force g is given by 



1 

z(Ag) 



E 7 =i h j (E,x 1 ,x 2 ) 



E 



£1,2 e 



-f3(E-gx 2 )_ 



-(3{E- gj x 2 ) ' 



(B-2) 



■/3(E-gx 2 )_ 



T,j = l hj(E,X!,X 2 ) 



where 
Z{(3,g) = J2 e ~ 

E,Xl ,X2 

(B-3) 

is an approximate partition function at <?, estimated us- 
ing Zfc, i.e., the partition functions at gk- The partition 
functions Zk are obtained self-consistently from the fol- 
lowing equation 



E 



,-/3(E-g k x 2 ) _ 



E,X\ ,X2 



Ej=i h J (E,x 1 ,x 2 ) 



(B-4) 



We iterate the above equation with starting values Zk = 
1 for all k. The convergence is monitored by estimating 
the amount of change after each iteration. We gauged it 
by calculating the quantity 



E 



7 {m) 



T {m-l) 



,(m) 



(B-5) 



where Z^ 1 is the value of Zk at the mth iteration 



When 

this quantity falls below some predefined quantity e , the 
convergence is achieved. We take e 2 = 1CU 14 . The force- 
distance isotherms obtained by this procedure in 1 + 1 
dimensions is shown in Fig. [2Kb). 
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